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LATTICE MODEL THEORY OF THE EQUATION OF STATE 
COVERING THE GAS, LIQUID, AND SOLID PHA ES 

ABSTRACT 

This studj' aims at obtaining the three stable state's of matter and the cor- 
responding phase tiansitions with a single model. Patterned after Lennard- 
Jones and Devonshires's theory, a simple cubic lattice model containing two fee 
sublattices (a and P) is adopted. The interatomic potential is taken to be the 
Ltnnard-Jones (6-12) potential. Employing the Cluster Variation Method, the 
Weiss and the pair approximations on the lattice gas failed to give the correct 
phase diagrams. Hybrid approximations are then devised to describe the lattice 
term in the free energy. A lattice vibration term corresponding to a free vol- 
ume correction is included semi-phenomenologically. The combinations of the 
lattice pa it and the free volume part yield the three states and the proper phase 
diagrams. To determine the coexistence regions, the equalities of the pressure 
and Gibbs free energy per molecule of the coexisting phases are utilized. The 
ordered branch of the free energy givei rise to the solid phase while the dis- 
ordered branch yields the gas and liquid phases. It is observed that the triple 
point and the critical point quantities, the phase diagrams and the coexistence 
regions plotted are in good agreement with the expeiimental values and graphs 
for Argon. It is also apparent in the investigation that the vibrational term was 
vital in getting the three phases. 



LATTICE MODEL THEORY OF THE EQUATION OF STATE 
COVERING THE GAS, LIQUID, AND SOLID PHASES 


I. INTRODUCTION 

At present, there is still a need for a unified theory that will cov#»** 'll three 
phases of matter as well as generate the proper transformations observed among 
phases. This study is an attempt to meet this need by utilizing a single model 
to formulate the equation of state of a simple substance semiquantitatively. 

The . add proposed here is the rigid lattice gas model of Horiguchi and Tanaka* 
with an effective vibration ol the atoms about the lattice sites to represent a 
free volume correction. Section II describes the model while Section III presents 
the results of the Weiss and the pair approximations in the Cluster Variation 
Method (CVM)^ as applied to the lattice part of the model. A new set of approxi- 
mations called the hybrid approximations is presented in Section IV, based on the 
CVM. The effective free volume correction is also formulated. Section V de- 
scribes the results and compares them with the results of the expandable model 
calculations by Mori, et al.^ and the experimental studies for Ai'gon.^ 


II. THE MODEL 

According to classical statistical mechanics, the Helmholtz free energy of a 
system of N particles is given by 

f = -V- “ 

where is the canonical p.iitition function of the system, is the Boltzma^m 
constant and T is the absolute temperature. The division of the free energy into 
the kinetic part 


K, = Nk^T - . 


with = (h^/ 27 , mk T)'’ ‘ , and the configurational part 


( 2 ) 
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where is the configurational integral 


Qn 


_L j e-^U(rN) 

N ! 


( 4 ) 


is the natural way of representing the total free energy. The evaluation of the 
integral in (4), however, is a formidable task. An alternative classification of 
the contributions to the free energy is offered by the lattice gas model which • 
avoids the difficult integrations. 

The total free energy is alternatively split into two quantities 

F ^ + F'' 

where F^ is the lattice part arising purely from the configurational free energy 
Fp and is evaluated conveniently in the Cluster Variation Method. The effective 
free volume F contains part of the configurational free energy and part of the 
kinetic free energy which are usually dropped in lattice model calculations. 

In this work, the rigid lattice model (RLM) is adopted because in the pressure 
range of interest, the lattice constant of the crystalline state does not change 
appreciably. Figure 1 shows a simple cublic lattice with f.c.c. a and /3 sub- 
lattices. N atoms are distributed over the L lattice sites forming the above 
lattice structure and with a fixed lattice constant a. The number density is then 
defined by 


P 



( 6 ) 


where v = V/Na'^ and V = La^. 

Similarly, the number densities of a sites and /3 sites are 


Pa 



(7) 


Pp=T- 


( 8 ) 
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where and Lo are the nunr.ber of the a sites and sites, respectively, vhile 
Na and are tne corresponding number of atoms on the a and sites. One 
has the following relations: 


Na + N;, - N. 

(9) 

Ljjj + = L, 

(10) 

, , 1 T 

(11) 

k - 9 L. 


A long range order (L.R.O.) para.neter is defined by 


R" (Pa -P)/(l -P). 


(12) 


Since the completely ordered state is described by a state in which all a sites 
are occupied, is equal to unity and hence K = 1 . In the completely disordered 
state, both the a and /3 sites are equally populated, hence, p^ = p^ = p and 
consequently R = 0. As defined by (12), therefore, the L.R.O. parameter is 
properly chosen. In terms of this parameter, the number densities on the a and 
/3 sites are written as 


P« = P + (1 - P) R 
P^ = P - (1 - P) R- 

It is of interest in this paper to investigate systems described by the Lennard- 
Jones potential 


0(r) 




(15) 


with a denoting the distance where the potential is zero and - e being the minimum 
value of the potential. The following dimensionless (reduced) quantities are 
introduced for calculational convenience. 


T = F^Z., g = i. (It* 

e e e € 

The free volume is now introduced as an effective vibration of the atoms about 
the lattice sites. 
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The free energy of a molecule confined in a harmonic potential well is known to be 


Tv = _X'enq=| + T'Cn(l 




( 17 ; 


where a? is a reduced frequency defined by 2; =tTo;/e ,oj being the classical fre- 
quency. The symbol q represents the canonical partition function of a simple 
harmonic oscillator. Assuming that all the molecules see identical confining 
cages (Einstein model), the vibrational free energy per lattice site may be written 
as 




(18) 


The quantity oj is an unknown function of the parameter R and/or density p. One 
knows by physical considerations that it should be an increasing function of these 
quantities since the stiffness of the confining potential increases with R and/or 
p . In this paper, a linear dependence on R is proposed, namely 

S = 2„K + a,, (19) 

where and oiq are adjustable parameters. The parameter can be inter- 
preted as a measure of the relative stiffness of the confining potential in the dis- 
ordered configuration compared to that of the ordered state. The parameter Zq 
on the other hand, can be chosen to give the correct triple point temperature. 
Equations (18) and (19) define the free volume term which will be added to the 
lattice free energy in later discussions. 


ni. THE WEISS AND THE PAIR APPROXIMATIONS 

In this section, the consequences of the purely lattice free energy are examined 
in the Weiss and the pair approximations of the Cluster Variation Method. 

The variational function for the one-site or Weiss approximation in the CVM 
is given by the trial free energy per lattice site. 


5 
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[30(1) + 60(2) + 40(3) + 30(4)] 


- [30(1) - 60(2) + 40(3) - 30(4)] (1 - p)^ r2 

+ iTl[p + (1 -p)R] in[p + (1 -p)R] + (1 -p) (1 -R)'en[(l - p) (1 - R)] 

+ [p - (1 -P) R] ln\p- (1 -p) R] + (1 -p) (1 + R)^n[(l - p) (1 + R)]} (20) 


in terms of the L.R.O. parameter R and U (k) =0 (k)/ e which is the reduced pair 
interaction with the kth nearest neighbor atom. Upon specifying the temperature 
^ anc . olume v, Eq. (20) becomes a function of R alone and minimization is 
easily done by setting 



(21) 


The non-zero solution for R yields the ordered state oi the free energy per 
lattice site^^j^p. The disordered state corresponding to the fluid state arises 
from Eq. (20) with the solution R = 0, i.e. 

^Dis = + ^^(2) + 4U(3) . 3U(4)] + T[p'&np + (1 + p) ^n(l - p)] . (22) 

The system's actual free energy is then deduced from 

(23) 

From thermodynamics, expressions for the pressure and Gibbs free energy 
per molecule (chemical potential) can be obtained. Specifically, for the rigid 
lattice model these are: 


P = 



(24) 


R = (T + P ) V. 


(25) 
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The ordered pressure and chemical potential follow from (24) and (25) if the 
ordered free energy is substituted while the corresponding disordered quantities 
are obtained if the disordered free energy is used. 

To avoid enormous graphical construction in determining the coexistence regions, 
the coexistence conditions for two phases A and B are used, namely 


P(P,) -- P(P^) 


(26) 




(27) 


for a given temperature T. This pair of conditions is equivalent to Maxwell's 
equal-area construction in the P-v isotherm or the common tangent construction 
in the? - v isotherm, where? = fxv denotes the reduced Helmholtz free energy 
per molecule. 


Calculations with the purely lattice term in the Weiss approximation as described 
above yield a couple of major difficulties: 

(1) It is not possible to get the liquid state and the liquid-gas critical point 
in the Lennard-Jones potential due to the stiff repulsion in the first 
nearest neighbor. A very soft core potential is necessary to obtain a 
fluid loop in the disordered branch of free energy. 

(2) Adopting an unrealistically soft core, a fluid loop is realized. However, 
the liquid appears in the wrong region of the phase diagram, namely 
gas, solid, liquid in that order with decreasing volume. 


The Weiss approximation, one may conclude, has favorable features it one wishes 
to describe remelting at high pressures where the existence of the vapor pressure 
curve and the critical point is not of concern. In the region near the triple point, 
however, which is the focus of this investigation, the Weiss approximation cannot 
give the proper phase diagrams even with the addition of a free volume correction. 

One order higher in the hierai’chy of approximations in the CVM is the two-site 
or pair approximation. The variational free energy of the lattice in this approxi- 
mation is given by 


7 
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?■- = 3<U(1) X^(l) + T.^[X,/1) in X^gO) 


* (/>a - in(p^ - X.,(l)) + (Pg - X.^(l)) in(Pg - X.^(l)) 


*0-Pa-Ps* X„/l)) in(l -P^-Pg* X./!))]) 


t 3{U(2)X.„(2) < Tl\<.(2)^u)^,,(2) 


t 2(/0„ - X.,(2)) in(p, - X,.(2)) + (1 - 2p. + X.„(2)) in(l - 2,=. + X„(2))) 


+ 3{U(2) Xa^(2) + T [Xg^(2) in Xgg(2) 


+ 2 (Pa - Xaa(2)) in(pg - Xgg(2)) + (l-2pgtXgg(2))in(l-2pgt Xgg(2)l } 

* 4{U(3) X.a(3) t T[)(,^(3) in X^g(3) 

* (p» - X.a(3)) in(p^ - X.a(3)) + (Pg - X.a(3)) in(pg - X„a'3)) 

+ (1 - Pa - Pp + x„a(3)) ^n(i -P^-Pg* x,a(3))l> 

+ I <U(4) X„(4) + T[X„„(4) in X„„(4) 

+ 2(Pa - X..(4)) in(j,^ - X„(4» + (1 - 2p„ t X,,(4)-tn(l - 2p, t X.„(4))]) 
+ 1 (U(4) Xaa<4) + TlXaa(4) in Xgg(A) 

*2(fig- Xgg(4))iv(pg - Xgg(4)) + (1 - 2pg + Xaa(4))-tn(l - 2pg + Xaa(4))) ) 

31 ~ 

- ~2 Pa ^ ~ Pa^> + + ^^>( 28 ) 


The notation (j) represents the pair correlation parameter betwe'm a pair of 
jth neighbors one of which is an a site, the other a yS site. The other pair param- 
eters are defined in the same manner. This trial free energy is minimized with 
respect to the six pair correlation parameters (1), X^^ (2), X^^ (2), X^^ (3), 
X^ (4) and X^^ (4), and with respect to R which enters through Pa and p^ . Tte 
minimum value with R non-zero gives rise to the ordered state free energy,^!, . 
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In analogy with the Weiss approximation, the disordered free energy in the 
two-site approximation is obtained with R = 0 and the pair parameters related 
byX..(k) = X3.(k)=X(k), 

f fejs = + T[X(1) 'enX(l) + 2{p - X(l)) U{p - X(l)) 

+ (1 - 2p + X(l))-en(l - 2p + X(l))]} 

+ 6{U(2) X(2) + T[X(2) In X(2) + 2(p - X(2)) ln{p - X(2)) 

+ (1 - 2p + X(2)) in{\ - 2p + X(2))]} 

+ 4{U(3) X(3) + T[X(3)^nX(3) + 2(p - X(3)) ^Tv(p - X(3)) 

+ (1 - 2p + X(3)) - 2p + X(3))]} 

+ 3{U(4) X(4) + T[X(4)^ia X(4) + 2(p-X(4)) ln{p - X(4)) 

' (1 -2p + X(4)) ln(l - 2p + X(4))]} 


- 3lT[p"fnp + (1 -p) -tnfl - p)] . 


(29) 


Tne pair parameters are determined by minimizing (29) with respect to the 
X(k), k = 1, 2, 3, 4. The ordered and disordered pressure and chemical potential 
are determined from (24) and (25) with the ordered and disordered free energies 
substituted, respectively. 

Calculations in the two-site cluster approximation showed a significant improve- 
ment in the sense that even with the Lennard-Jones potential, the fluid loop can 
exist in a reasonable range of temperature. This implies that if one disregards 
the presence of the solid phase in the diagram, a finite critical temperature is 
realized, making the existence of a liquid state more feasible. 

A basic difficulty however is the instability cf the results with variation in the 
order of the nearest neighbors included in the calculation. If only the nearest 
and second neighbor pairs are included, the free energy of the ordered branch 
is lower than the disordered branch of free energy. However, when third or 
fourth nearest neighbors are included, the disordered branch becomes lower. 
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It is in the spirit of these difficulties in the Weiss and pair approximations that 
a physically reasonable way of combining them was devised as described in the 
following section. 


IV. THE HYBRID APPROXIMATIONS 

The hybridization procetlure is based on a sound physical observ.’tion that 
neighbors which are situated farther from the site under consideration are more 
blurred in their details than the nearer sites. Applying this to the lattice problem, 
the first nearest neighbors are treated in the two-body approximation while the 
farther neighbors (2nd, 3rd, 4th, etc.) are described in the Weiss approximation. 
Variation? of this process are clearly possible. 

It is convenient in the hybrid approximations to adopt the notation [2 - (m)l] to 
indicate that the effect of the nearest neighbor interaction is calculated in the 
pair approximation while the next m shells of interacting neighbors are treated 
in the Weiss approximation. Hybrid approximations [2-l] , [2-(3)ll, [2-(5)l], 
and [2-(7)l] are investigated. The trial free eneigies in these hybrid approxi- 
mations are the following: 

?»'[2 - 1] = 3[U(1) +U(2) (pj ^ 3T[X,^(l)'tnX,^(l) + (p, 


+ X,^(1))J -p{p^ Inp^ + (1 -pj + (1 -p^)^n(l (30) 


T 42 - (3) 1] = fL[2 - 1] . 4U(3)p^p^ + |U(4) (p.J + pj) 

T42 - (5) 1] =T -[2 - (3) 1] + 12 U(5)p^p^ + 6U(6) (pj + p2) 
[2 - (7) 1] = ?42 - (5) 1] + 3U(7) (p2 + p2) + 4U(8) p^p^. 


(31) 

(32) 

(33) 


These trial free energy expressions are minimized with respect to R and the 
pair function (1) giving rise to the ordered values of the free energies in the 
different h}'brid approximations for non-zero R, as functions of density P and 
temperature T. 
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In order to obtain the disordered branch of the free energy in these approxima- 
tion, R is set to zero and the resulting expressions are minimized with respect 
to the single pair parameter X(l). One then gets the following: 

TLjg[2-l] = 3 [U(l)X(l) + 2U(2)p2] . 3T[X(l)^nX(l) + 2(p-X(l))H:n(p - X(l)) 


f (l-2p+X(l)'en(l -2p + X(l))-5T[p^np + (l-p)'fcn(l -p)] (34) 

[2- (3) 1] = rfeis [2 - 1] + [4U(3) + 3U(4)] ^35) 

*^Dis - (5) 1] = “ (3) 1] + [12U(5) + 1 2U(6)] p2 ' (36) 

^Dis ^2 - (7) 1] = t2 - (S) 1] + [6U(7) + 4U(8)] p^. (37) 


The corresponding ordered and disordered pressure and chemical potential follow 
from (24) and (25). 

Ail the above hybrid approximations on the lattice Helmholtz free energy per site 
displayed basically the same characteristics when plotted as functions of volume. 
A finite value for the critical temperature T^. ^s a common feature reminiscent 
of the pair approximation. Although the lattice part alone does not suffice to 
correctly describe the thermodynamic phase diagrams and the three states, all 
the hybrid approximations resulted into a lower ordered free energy consistently. 
They are therefore susceptible to a corrective treatment which essentially raises 
the ordered relative to the disordered free energy like the free volume correction 
described earlier. 

There are four combinations arising from the four hybrid approximations on the 
lattice investiga-.ed. The trial free energies of the combined lattice and free 
volume parts are as follows: 

T[2 - 1] =?L12 - 1] ♦ f'' (38) 

? [2 - (m) 1] = [2 - (m) 1] + T'' m = 3, 5, 7. (39) 

where7^ is defined by (18) and (19). 

As before, the ordered state free energy in the various approximations are 
obtained by minimizing the above equations with respect to R and X^(l). The 
corresponding disordered free energy branches follow by setting R = 0 and 
minimizing the corresponding exp**essions 
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T [2 - 1] = [2 - 1] + (40) 

r „,5 12 - (m) I) = [2 - (m) 1] + r»,j m = 3. 5. 7 (41) 

with respect to X(l). Use of equations (24) and (25) enables one to get the expres- 
sions for the ordered pressure and chemical potentials 


P [2 - ll = - f 

ORD ^ J ORD 


, /3^o«d [2 - 11\ 


(42) 




3^okd[2-(">)11 




®ORD 


ORD 


ORD 


T.L 

m = 3, 

5. 7 

(43) 

[2 - 

1]) v 


(44) 

) V 

tn = 3, 

5, 7 

(45) 


and the corresponding disordered pressure and chemical potentials. 


In all approximations, the system's actual free energy is given by the lower of 
the two branches of free energy, namely 


^“Min{fojjp, fpis}- 

The results uf these approximations are presented and discussed in the next 
section. 


(46) 


V. RESULTS AND CONCLUSIONS 

Having obtained the free energy and pressure expressions in the hybrid approxi- 
mations with free volume correction, one makes use of the coexistence conditions 
(26) and (27) to obtain the phase transformations and the coexistence regions for 
melting, ondensation and sublimation. Applying these conditions appropriately 
over a range of temperature T, the PT diagrams shown in Figure 2 for a simple, 
classical system are generated for various hybrid approximations with free 
volume correction. The coexistance diagram is shown in Figure 3 for the I 2 - (7)l] 
hybrid approximation with lattice vibration. The vibrational parameter is taken 
to be .01 for all hybrid approximations while zo = 1.256 for the [2 - (3)1] and 
[2 - (5)1] approximations and z^ = 1.415 for the [2 - (7)1] h 3 rbrid approximation. 


12 



14 


The typical behavior of an isotherm for temperatures above the triple point is 
shown in Figures 5a, b, and c where an isotherrr of temperature T = 0.71 is 
plotted in the pressure, Helmholtz free energy per molecule, and chemical po- 
tential against density diagrams in the [2 - (5)1] approximation with free volume 
correction. 

The typical behavior of an isotherm for temperatures below the triple point tem- 
perature is shown in Figures 6a, b, and c. Plotted in the pressure, Helmholtz 
free energy per molecule, and chsmical potential against diagrams is the T = 0.68 
isotherm in the [2 -(3)1] approximation with free volume correction. There is 
a slight variation here since two tr.''nsitions, gas to liquid and gas to solid, com- 
pete because both are possible under the criteria (26) and (27). It is the gas to 
solid transition which is realized, however, since it corresponds to a lower 
Helmholtz free energy per molecule and lower chemical potential. 

In Table 1, the triple point quantities tre shown. The results of the present 
rigid lattice calculations (Row I) are compared with the expandable lattice calcu- 
lations (Row n) by Mori et al.^ and the experimental values for Argon (Row IH). 
Similarly, the critical point quantities are displayed in Table 2. Within the 
Jiermodynamic region of interest, it is apparent that the present rigid lattice 
’heory is better than the expandable lattice theory. Calculations for the latent 
heat of transformations using the Clausius -Clapeyron equation and the phase 
diagrams obtained in the present theory are being done. The authors hope to 
present the results in the future. 

To improve the present theory, it is desirable to obtain a more accurate 
formulation of the free volume term, ultimately avoiding the phenomenological 
part. Raising the order of the CVM approximation is a possible option. A better 
alternative way is suggested by introducing more types of lattice sites, effectively 
shrinking the lattice spacing. In this process, the effect of vibration is gradually 
and partly incorporated into the lattice term. A physical limit in the spacing 
might exist after which the ad hoc addition of the vibrational or free volume 
term becomes unnecessary. 

The nature of the convergence of the results in the hierarchy of approximations 
in the CVM toward the correct ones is not yet clear (even for systems related 
to the lattice gas, e.g. magnetic systems, binary alloys, etc.). A related question 
is the stability of the general results with respect to the order of neighbors 
included. As the interaction becomes long-ranged, one must naturally include 
more neighbors in the calculations. The authors believe that the order of the 
approximation and the range of interaction are intimately connected. The hybrid 
approximations introduced in this paper are the simplest ways of investigating 
the effects of longer-ranged interactions. It is obvious that higher orders in the 
CVM can also be hybridized v/hich, hopefully, may shed more light to the present 
querry. 
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Q ^ SITE 

Figure 1. The simple cubic lattice formed by o- and /3 sites. Both the ^ and 
/3 sublattices are face-centered cubic in structure. 
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<J.t 0.7 0.? 0.1 /O 1.1 i.i / 

T 

Figure 2. The pressure vs. temperature (PT) diagrams in the hybrid 
approximations with free volume correction (wq = .01). 





Figure 3. The coexistence diagram in the pressure vs. density (P - p) 
plot and in tne [2-(7)Il approximation with free volume correction 
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Figure 4. The liquid-gas coexistence curves. Shown here are the [2-1] , 
[2-(5)l],and [2-(7)l] approximations, and the experimental data ’joints 
for Argon. Since the liquid and the gaseous phases are both derived 
from the disordered branch of free energy, the vibrational parameters 
are unimportant in this diagram. 
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Figure 5b. The isotherm T = 0.71 in the [2 - (5)1] approx- 
imation with free volume correction (u)q = .01, Zq = 1.256); 
Helmholtz free energy per molecule vs. density (T„ ->o) 
diagram. 
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Figure 5c. The ieotherm T = 0.71 In the [2 - (5)1] approx- 
imation witli free volume correction (u)q = .01, Zq = 1.256); 
chemical potential vs. density (g-g) diagram. 
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Figure 6a. The isotherm T = 0.68 in the [2 - (3)l] approx- 
imation with free volume correction (ujq = ,01, = 1.256) 

plotted in the pressure vs. density (P - P) diagram. 
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T » 0.68 

[2-(3)lj APPROXIMATION 
ORDERED 



Figure 6b. The isotherm T = 0.68 in the [2 - (3)lj approx- 
imation viui free volume correction (u),) = .01, Zq = 1.256) 
plotted i^thfc Helmholtz free energy per molecule vs. 
density (f^ - p) diagram. 
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Figure 6c. The isotherm T = 0.68 in the [2 - (3)l] approx- 
imation with free volume correction (ujq = .01, = 1.256) 

plotted in the chemical potential vs. density (g - c) diagram. 
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Table 1 

The Triple Point 


Model and 
Approximation 

T. 

r 

P* = 2P^ 

m 

P, V, /T, 

R 

(LRO) 

1. Rigid 




■■■I 



(a) [2 -(3)1] 

(G) 




0.8364 

0 

Zq = 1.256 

(L) 

0.713 

0.0198 


0.0384 

0 

(Vq = .01 

(S) 



1.0010 

0.0278 

0.9981 

(b) [2 -(5)1] 

(C) 



159.52 

0.9571 

0 

Zo - 1.256 

(L) 

0.700 

0.0042 

1.187 

0.00712 

0 

= .01 

(S) 



1.00024 

0.00600 

0.99960 

(c) [2 -(7)1] 

(G) 




0.9742 

0 

z. = 1.415 

(L) 

0.698 

0.0032 


0.00536 

0 

o 

II 

(S) 



1.00020 

0.00458 

0.99964 

n . Expandable 







(a) Without 

(G) 




0.7686 

0 

Free Vol. 

(L) 

1.487 

0.055 


0.0676 

0 

MC2(1) 

(S) 




0.0463 

0.9997 

(b) With Free 







Volume 







ME4 (1 , 1 ) 







1 I 

(G) 




9.9206 

0 

(2.10) 

(L) 

0.7789 

0.00550 


0.0101 

0 


(S) 




0.0079 

1 


(G) 



26.57 

0.7320 

0 

(3,12) 

(L) 

0 8509 

0.02344 

1.654 

0.0456 

0 


(S) 



1.170 

0.0323 

1 

III. Experiment 

(L) 

0 700 

0.00151 

1.186 

0.00256 


(Argon) 

(S) 



1.035 

0.00223 
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Table 2 

The Critical Point 


Model and 
Approximation 

T 

C 

P* = 2P 

C C 

m 

P V /T 

c c c 

T./T, 

I. Rigid Lattice 
(a) [2-1] 


0.074 

3.20 

11 

II 

[b) [2-(3)l] 


0.1080 

3.16 

■m 


(c) [2-(3)ll 

1.237 

0.1382 

3.25 

0.363 


(d) l2-(7)ll 

1.282 

0.1431 

3.26 

0.364 

1.837 

n. Expandable 

(a) Without Free 
Free 
Volume 

MC2(1) 

2.049 

0.2001 

3.682 

0.3596 

1.378 

MC2(a) 

1.751 

0.1271 

4.023 

0.2921 


(b) With Free 
Volume 

ME4 

=0 

2.2974 

0.3025 

2.665 

0.3509 


=2 

1.5331 

0.2238 

3.394 

0.4956 

1.968 

CO 

II 

V 

1.2025 

0.1999 

2.920 

0.4852 


in. Experiment 
(Argon) 
e/k = 119.8“ K 

1.26 

0.116 

3.16 

0.291 

1.80 

cr = 3.405 A 
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